RATIONAL APPROXIMATIONS FOR THE QUOTIENT OF 

GAMMA VALUES 



KH. HESSAMI PILEHROOD 1 AND T. HESSAMI PILEHROOD 2 

Abstract. In this paper, we continue to study properties of rational approximations 
to Euler's constant and values of the Gamma function defined by linear recurrences, 
which were found recently by A. I. Aptekarev and T. Rivoal. Using multiple Jacobi- 
Laguerre orthogonal polynomials we present rational approximations to the quotient 
of values of the Gamma function at rational points. As a limit case of our result, 
we obtain new explicit formulas for numerators and denominators of the Aptekarev 
approximants to Euler's constant. 



1. Introduction 

Recently, A. I. Aptekarev and his collaborators in a series of papers [U [TT] proposed 
a sequence of rational approximations p n /q n converging to Euler's constant 

7= lim (1 + - + -•• + - -In A; ] . 

fc^oo \ Z k J 

sub-exponentially. The numerators p n and the denominators q n of the approximations 
are positive integers generated by the following recurrence relation: 

(16n - 15)g n+ i = (128n 3 + 40n 2 - 82n - 45)g n 
1 1 ' - n 2 (256n 3 - 240ra 2 + 64n - 7)q n ^ + n 2 {n - l) 2 (16n + l)g n _ 2 

with initial conditions 

/ 2 \ Po = °> Pi = 2 ' V2 = 31, 

1 ' Qo = 1, qi = 3, q 2 = 50 

and having the following asymptotics: 



e V2n 



«» = < 2 "> ! ^(vw +0( ' rl/u 



(3) Pn ~ ldn = (2n)!-^- + 0(n-W] 
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In [BJ, explicit formulas for the sequences p n and q n were given in terms of hypergeo- 
metric sums 



n / \ 2 

n 



(4) fc = EL ( n + k v-> 

k=0 ^ ' 



where i7 n is an n-th harmonic number defined by 

1 

k + a 



Hn{a) = ir T ^—, n>l, H (a) = 0, H n = H n (0). 

< * ic -4- rv 



k=l 



These representations provide a more profound insight on arithmetical properties of 
these sequences that can not be deduced from the recurrence equation (Op). From (Tjfl) it 
follows that q n and p n are integers divisible by n\ and n\/D n respectively, where D n is 
the least common multiple of the numbers 1,2, ... ,n. This implies that the coefficients 
of the linear form fl3]) can be canceled out by the big common factor jj-, but it is still 
not sufficient to prove the irrationality of 7. The remainder of the above approximation 
is given by the integral [3] 



00 

g(0,0,D (x) e - ln(x) dx = pn _ lgn 

where for a%, a 2 > — 1, b > 0, 

1 2 d n 

(5) Q^ b \x) = - J] W-rJOiA (1 - x 



2 

2n 

n! 2 - LJ - J c?x r 



is a multiple Jacobi-Laguerre orthogonal polynomial on Ai = [0, 1] and A 2 = [1, +00) 
with respect to the two weight functions 

Wl ( x ) = x ^{\ - x )e- bx , w 2 (x) = x a2 (l - x)e~ bx 

satisfying the following orthogonality conditions [2]: 

(6) / x"Q ( ^ 1 ' a2 ' b) {x)w k {x)dx = 0, j,k = 1,2, v = 0,l,...,n- 1. 



Note that, as the values of the parameters a\, a 2 tend to and 6=1, the limit values of 
Qn'°' 1 \x) satisfy the orthogonality relations (J6j with the weights W\{x) = (1 — x)e~ x , 
Wi{x) — (1 — x) ln(x)e _x . 

In [6J, the authors showed that for a. b <G Q. a > — 1, 6 > the integrals 

L2n+a+l ^00 

(7) - ■ - / Q^ b \x)x a ln(x)e- bx dx = p n (a, b) - (ln(6) - V(a + l))g„(a, b), 



T(a + 1) 

where if)(z) = T'(z)/T(z) defines the logarithmic derivative of the gamma function, pro- 
vide good rational approximations for ln(6)— ip(a+l), which converge sub-exponentially 



2 



to this number 

Pn{a, b) _ _ ^ + = 2ne -2V2bn^ + o^- 1 /?)) as n oo, 



9n(a,6) 

and gave explicit formulas for the denominators and numerators of the approximations 

2 



n / \ 2 

g n (a,6) = ^rj (a + l) n+fc 6 n - fc , 

k=0 ^ ' 
n / \ 2 

(9) P«M) = X)(l) ( a + 1 W b n ~ k H n+k {a) 

1,—n \ / 



fc=0 

n A; n— fc 



\k + I J \m J \ I 

k=l m=0 1=0 v / \ / \ 



here (a) n is the Pochhammer symbol defined by (a)o = 1 and (a) n = a(a + 1) • • ■ (a + 
n — 1), n > 1. 

Recently, T. Rivoal [9] found another example of rational approximations to the 
number 7 + ln(x), x € Q, x > 0, by using multiple Laguerre polynomials 



A n (x) = ?-(x n {x n e-*)^)^ = e x - 2 F 2 ( n + n + 1 
nr \ 1,1 



.2; 



(for the validity of the last equality see [7J §4.3]). Here p F g is the generalized hyperge- 
ometric function defined by the series 



F 

o ± q 



Oil , . . . , Olp 



(ai)fc ■ • • (a p )fc z k 

(A)* • • • (A,)* W 



E 



which converges for all finite 2 if p < q, converges for \z\ < 1 if p — q + 1, and diverges 
for all 2, 2 7^ if p > q + 1. The series p F q terminates and, therefore, is a polynomial in 
z if a numerator parameter is a negative integer or zero, provided that no denominator 
parameter is a negative integer or zero. 

The polynomial A n (x) as well as the hypergeometric function 2-^2 ( ri + 1 ' ri + 1 1 — x) sat- 
isfies a certain third-order linear recurrence relation with polynomial coefficients in n 
and x. The polynomials A n (x) are orthogonal with respect to the weight functions 
tO\{x) = e~ x , U2(x) = e~ x ln(a;) on [0, +00), so that 

POO 

(10) / x v A n {x)uok{x) dx = 0, A; = 1,2, v = 0, 1, . . . , n — 1. 

This implies that A n (x) is a common denominator in type II Hermite-Pade approxima- 
tion for the two Stieltjes functions 

e X {x) = r — dt, s 2 (x) = r e ~ l in(t) dt 

J x t Jo x t 
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at x = oo. The remainders of this approximation are given by the integrals 



Ri,n( x ) :z 



x — t 



e~ l dt = A n {x)£ 1 {x) + B n (x) 



A n (t) 
x — t 



ln(t)e~* dt = A n (x)£ 2 (x) + C n (x), 



which in view of the orthogonality relations (1101) satisfy the same recurrence as the 
polynomial A n (x) (see [HI §5.2]). By linearity, it is easily seen that the remaining 
polynomials B n (x),C n (x), and the function 



R n (x) = \n(x)Ri <n (x) - R 2 , n {x) 



. , . ln(x) — m(t) , , 
x — t 



are also solutions of the same recurrence. Suppose that this third-order recurrence has 
the form 

U n +3 — S n U n+ 2 + tnUn+l + r n U n . 

Then rational approximations to ln(x) +7 are given [9] by a so called "determinantal" 
sequence 

A n (x) R n {x) 



11) SJx) 



A n+l (x) R n +i{x) 



Q n (x)(\n(x)+^)-P n (x), Q n (x),P n (x) G Q[x] 



which along with the sequences Q n (x) and P n (x) satisfies a "determinantal" recurrence 
of the form 

The sequence P n (x)/Q n (x) converges to \n(x) + 7 sub-exponentially 
Pjx) 



ln(x) + 7 



Qn(x) 



< Co ( x ) e - 9 / 2 - 1/3 ™ 2/3 + 3 / 2 - 2/3 " 1/3 , Q n { x ) = (e 3 * 1/3n2/3 -* 2/3nl/3 ), 



and for denominators of coefficients of the polynomials P n and Q n one has 

n\ Q n (x) G Z[x], D n n\ P n (x) G Z[x] 

(the first inclusion was formulated in [9] as a denominator conjecture, and both inclu- 
sions were proved in [6]). 

In [10], the similar scheme was applied to the Laguerre polynomials 



A nA (x) = ■^(x n - a (x n+ °e- x )W)W 



(a + 1) 



e x 2 F 2 



nl 



rv. 



a + n + 1, n + 1 
a + 1,1 



which are orthogonal on [0, +00) with respect to the weights e~ x and x a e~ x , to get 
rational approximations for T(a + l)/x a , where a, x G Q, a > —l,x > 0. In this case, 
rational approximations are defined by the "determinantal" sequence 



(12) S n (a,x 



Q n (a, x) T(a + 1) — P n (a, x) x a , 



where 



Rn.ai^^} x R 



A n+ i A (x) Rn +1) a(x) 

Rn,a,a\ x ) ari d Rn,a,c\ x ) = 

JO 



n,a,0\ x ) 



X — t 



The sequences P n (a, x), Q n (a, x),S n (a, x) are solutions of the third-order recurrence 

C 3 (n, a, x)V n+3 + C 2 {n, a, x)V n+2 + Ci(n, a, x)V n+1 + C (n, a, z)V n = 0, 

where the coefficients Cj(n, a, x), j = 0, 1, 2, 3 are polynomials with integer coefficients 
in n,a,x of degree 16 in n. Moreover, one has [TD] 

\Q n (a, x)T(a + 1) - P n (a, x)x a \ < Cl (a, x) e-^^n^+i/^^ya^/s^ 

Q n (a,x) = (e 3a;1/3 " 2/3 -^ /3ril/3 n 2a / 3 - 2 ) 
with the following inclusions proved in [6]: 

(n + l) 2 • n! /if +1 (den x) n+1 P n (a, x), (n + l) 2 ■ n\ /4" +1 (denx) n+1 Q n (a, x) G Z, 
where 

(13) Ai» = (dflno) B . J] 

p|dena 

and den a e N is the denominator of a simplified fraction a. 

In this paper, we continue to study properties of the above constructions and get 
rational approximations for the quotient of values of the Gamma function at rational 
points. The main result of this paper is as follows. 

Theorem 1. Let a\, a 2 , b e Q, a% — a 2 ^ Z, a%, a 2 > — 1, b > 0. For n = 0, 1, 2, ... , 

define a sequence of rational numbers 

q n (ai,a 2 ,b) = 2^y fc n _ fc J(a 2 + l) n+fc & , 

T/ien 

Va 2 - ai V 2 a n 2 (ln{ai,a 2 ,b) e n\Z[b], 
where /x™ defined in (TT3"j) , and t/ie following asymptotic formulas are valid: 

(14) g " (ai ' ° 2 ' 6) r(ai + !)/^ - 9 " K fll ' h) = ^ ! nVM^)/2 ( Cl + °^ 1/2 ^' 

g n (ai, a 2 , 6) = (2^)! ^ 1/4 _ (ai+aa)/2 (c 2 + 0(n~ 1/2 )) asn^oo 

with 

n . 2 sin(7r(a 2 - a0)r(a 2 + 1) 2 ^ +a ^/ 2 - 3 / 4 



&a 2 -«ir( ai + 1) Z ' ' 6 1 /4+(a 1 -a 2 )/ 2e 36/8 v ^ r ( a2 + 1 )- 

The linear forms ( IX4jl do not allow one to prove the irrationality of the quotient 
T(a 2 + l)/(r(ai + l)b a2 ~ ai ), however they present good rational approximations to this 
number 

T(a 2 + l)/b a * _ g w (a 2 ,ai,6) = ci g _ 2v ^ ^ 

For an overview of known results on arithmetical nature of Gamma values, see [10] and 
the references given there. 



Application of Zeilberger's algorithm of creative telescoping [SJ Ch. 6] shows that 
both sequences q n (ai, a 2 , b) and g n (a 2 , ai, b) satisfy a third-order linear recurrence with 
polynomial coefficients in n,ai,a 2 ,b of degree not exceeding 10 in n. We do not write 
down this recurrence equation for arbitrary a±, a 2 , b, since it is quite complicated in its 
general form. We indicate this recurrence only in several particular cases. 

Setting b = 1, a 2 = —1/2, ai = —2/3 and then a\ = —3/4 we deduce the following 
results. 

Corollary 1. The sequences 

*-£(-:*) - -eC ;*)(:: -Ml 

— /c — 

satisfy the third- order linear recurrence 

7776(n + 3)(n + 2)(1728n 4 + 6228n 3 + 7986n 2 + 4538n + 963)/„, +3 
= 36(2985984n 7 + 33654528n 6 + 155405952™ 5 + 379069848n 4 + 525871008™ 3 
+ 415828470n 2 + 175094660n + 30377295) (n + 2)/ n+2 

- 6(6n + ll)(6n + 13)(995328n 8 + 11384064™ 7 + 55911168™ 6 + 154077840n 5 
+ 260388504n 4 + 275942444n 3 + 178728654n 2 + 64534088n + 9914615) / n+ i 
+ (6n + ll)(6n + 5)(2n + l)(3n + l)(6n + 13)(6n + 7)(1728n 4 + 13140n 3 
+ 37038n 2 + 46106n + 21443)/ n 

with the initial values po — 1, pi = 43/54, p 2 = 12871/1458 and go = 1, Qi — 29/24, 
g 2 = 5149/384, and provide rational approximations for the following quotient: 

r(l/2) p» 2 W +Q(n -i/a^ as n^oo 

r(i/3) g n 6 lr(i/3) +ul V s 

Corollary 2. T/ie sequences 



fc=0 v 7 v 7 fc=o v 7 v 7 

satisfy the third- order linear recurrence 

2048(n + 2)(n + l)(2048n 4 - 1024n 3 - 304n 2 + 416n - 117)/ n+2 
= 32(4n + 5)(n + l)(262144n 6 + 753664n 5 + 331776n 4 - 365824n 3 + 17168n 2 
+ 76032n - 32229)/ n+ i - 2(4n + 3)(4n + 5)(4n + l)(524288n 7 + 1572864n 6 
+ 1486848™ 5 + 389120n 4 - 245824n 3 - 117296n 2 - 604n - 537)/ n 
+ (An + 3)(An - l)(2n - l)(An + 5)(4n + l)(4n - 3)(2048n 4 + 7168n 3 
+ 8912n 2 + 4928n + 1019)/„_! 



with the initial values p — 1, p\ — 37/64, p 2 = 50685/8192 and q — 1, qi — 19/16, 
g 2 = 6525/512, and provide rational approximations for the following quotient: 

r(i/4) g „ vr(i/4) +tA v 5 °°- 

Setting a,\ = 0, a 2 = a, b = 1 in Theorem [T] we get rational approximations for values 
of the Gamma function. 

Corollary 3. Let a e Q \ Z, a > —1. T/ien i/ie sequences 

^=^t a )^ n z a 1 ) {n+k) - and ^=g( n r)(n-^) (a+iw 

satisfy the third- order recurrence 

A 2 (n, a)/n+2 + Ai(n, a)/ n +i + A (n, a)/ n + A_ x (n, a)/„_i = 0, 

w/iere 

A 2 (n, a) = (n + 2)(16n 3 + 20n 2 a + 10na 2 + 2a 3 + 17n 2 + Una + 3a 2 + n - a), 
X^n, a) = -128n 6 - 224n 5 a - 112n 4 a 2 + 20n 3 a 3 + 42n 2 a 4 + 16na 5 + 2a 6 - 808n 5 

- 1228n 4 a - 601n 3 a 2 - 45n 2 a 3 + 47na 4 + 11a 5 - 1910n 4 - 2446n 3 a 

- 1085n 2 a 2 - 176na 3 - 3a 4 - 2035n 3 - 2063n 2 a - 715na 2 - 91a 3 - 887n 2 



592na - 113a 2 - 82n + Ua 



A (n, a) = (n - a + l)(n + a + l)(256n b + 576n 5 a + 544n 4 a 2 + 272n 3 a 3 + 72nV 

+ 8na 5 + 1296n 5 + 2384n 4 a + 1724n 3 a 2 + 596n 2 a 3 + 92na 4 + 4a 5 + 2448n 4 
+ 3518n 3 a + 1840n 2 a 2 + 404na 3 + 30a 4 + 2185n 3 + 2333n 2 a + 817na 2 + 93a 3 
+ 923n 2 + 668na + 125a 2 + 146n + 58a), 
A_i(n, a) = -(n - a)(n - a + l)(n + a + l)(16n 3 + 20n 2 a + lOraa 2 + 2a 3 + 65n 2 
+ hAna + 13a 2 + 83n + 33a + 34) (n + a) 2 , 

with the initial conditions 

Po — 1, p\ = 3 + a, p2 = 50 + 33a + 7a 2 , 

g = l, qi = ( a + l)(3-a 2 ), g 2 = (a + l)(a + 2)(a 4 + 2a 3 - 12a 2 - 11a + 50)/2, 

and provide rational approximations for the following value of the Gamma function: 

r( a + l)-^ = e - 2 ^f , 2 ™ +0(n- 1 / 2 )) as n -> oo. 
g„ Vr(l-a) / 

Another positive feature of our considerations is the fact that we can prove a simpler 
formula for the numerator p n (a, b) of rational approximations (jHJ) than a representation 
given by OH])- 



Theorem 2. For the sequence p n (a, b), n = 0, 1, 2, . . . , defined by Qj, a more compact 
formula is valid 

p n (a, b) = ( J (a + l) n+ k b n ~ k (H n+k (a) + 2H n _ k - 2# Jfc ). 

As a consequence, we get a simple formula for the numerator p n of rational approxi- 
mations ([3]) to Euler's constant. 

Corollary 4. TTie sequences p n and q n , n = 0,1,2,..., defined by the recurrence ([T]) 
with initial conditions (T2J) are given by the formulas 

n / \ 2 

fc=0 ^ ' 
n / \ 2 

P» = Z)(j {n + k)\(H n+k + 2H n . k -2H k ). 



2. Analytical construction 

Let ai,a 2 ,b be arbitrary rational numbers such that a± — a 2 Z, ai, a2 > —1, b > 0. 
The construction of rational approximations to the numbers r(a 2 + 1)/ (r(ai + l)6 a2 ~ ai ) 
is based on the integral 

poo 

(16) R n := R n {a u a 2 , b) = b 2n+1 / (x a2 - x ai )Q^' a2 > b) (x) e - bx dx, 

Jo 

where the polynomial Qn 1,a2,h \x) is defined in ([5]). 

Lemma 1. Let ai, a 2 , b G Q, a\ — a 2 Z, Oi, a 2 > — 1> & > 0. T/ten we /icwe 
i? n (ai, a 2 , 6) = q n {a u a 2 , b) — q n (a 2 , a 1 , by 



where q n (ai,a 2 ,b) is defined in Theorem^ 

i/2n+a 2 + l roo 

(17) q n ( ai , a 2 , b) = - / x a2 Q^ a2 ^e~ bx dx 

T(a 2 + 1) J 

and Ma 2 -aX" 9n(ai, a 2 , 6) G n! Z[&]. 

Proof. Substituting (jSJ) into (1TB]) we get 



18 

^2n+l roo _^ai-a2 



R n {a 1 ,a 2 ,b) = — — / dx 

n\ In 1 x 



n\ 2 J n 1 — x 

8 



(x n+a2 - ai (x n+ai (l - x ) 2 «+ 1 e -fe)W)(«) 



Denoting the first integral on the right-hand side of ( TTBj) by Ii and applying n-integrations 
by parts, we have 

( ^\njj2n+l roc ^n+a^—ax 



h 



IV. 



■f -{x n+a Hl - xf n+l e- bx ) {n) dx 

(1 - x) n+1 v v ' ' 



u2n+l poo / n+a 2 -a 1 \ ( n ) 

' ' —) x n+a ^l-x) 2n+1 e- bx dx. 

\n+l / ^ ' 



rv. 



(l-x) 



Now using the following relation for hypergeometric functions (see [H §3.4, (11)]): 
(19) 



dx 1 



\x a {l-xf] = (a-n + l) n x a - n (l-x) l3 - n 2F 1 



-n, a + 1 + f3 — n 
a + 1 — n 



x , n - a 4 N, 



with a = n + a 2 — cti, (3 = —n — 1, we get 



x 



n+a,2—ai \ ( n ) 



(1 - x) n+1 
and therefore, 



(a a -ai + l) n x Q2 - ai (l - x)- ln ~\F 1 



-n, a,2 — en — n 
a-2 — ai + 1 



X 



n\ Jo \ 



—rij <i2 — a\ — n 
«2 — «i + 1 



x 1 e 6:r dx 



(20) 



6 2n+1 (l + a 2 - ai) n y-v (-n) fe (a 2 - ai - n), 



rv. 



k=0 



k\ (1 + a 2 - ai) fc 



.r 



- / 
6 2n+l ^ ( 

k=0 ^ 



n + a\ — a 2 \ f n + a 2 — aA r(n + a 2 + + 1) 



fc 



n — k 



Thus we have 



r(a 2 + 1) 

6 a 2 



which implies ( ITTj) . Similarly, denoting the second integral on the right-hand side of 
f TTgj) by J 2 and integrating by parts, we have 



(21) 



(_1) n 5 2n +! Z" 00 / x ai ~ a2 \ 



(n) 



Applying formula (fl9l with a = a% — a 2 , /3 = — 1 we get 



x 



a±—a2 \ ( n ) 



1 — X 



( ai -a 2 -n + l) n x ai - a2 ~ n (l - xy n - L 2 F! 



-n, a\ — a 2 — n 



ai — a 2 — n + 1 

Substituting the n-th derivative into the integral (121 p and integrating by parts we obtain 
_ b 2n+l {a 1 -a 2 -n + l) n 

i5 — 771 



n 



1 2 



(1 - x)-™- 1 ,^ 



— n, ai — a.2 — n 
Q-i — 0-2 — n + 1 

9 



(«) 



x n+ai (l- x) 2n+l e- bx dx 



Using the following relation (see [TJ §3.4, (10)]): 



dx n 



;i - x) a ^~ s 2 F 1 



foe, P> 


x^j 


V * 





(5 - a) n (6 - ft), 

(*)n 



'l-x) a+p - 5 - n 2 F l 



a, (3 
5 + n 



x 



with a = —n, f3 = a\ — a 2 — n, 5 = a x — a 2 — n + 1 we can compute the n-th derivative 
in the last integral 



—n, a% — a 2 — n 
a\ — a 2 — n + 1 



x 



(») 



X 2*1 



(oi - a 2 + l) n n! 
- a 2 - n + l) n 
— n, a\ — a 2 — n 
ai — a 2 + 1 



-2n-l 



to get 



^(ai-os + l),, 



n! 



,r 



n+ai 



-n, ai — a 2 — n 
a 1 — a 2 + l 



x e 



x 



- bx 



dx. 



Comparing this integral with f )20|) we obtain 



r(ai + l) 

b ai 



and the desired expansion of R n (ai, a 2 , b) is proved. 

From the definition of the sequence q n { a ii a 2, b) it follows that q n (ai, a 2 ,b) G Q[b). To 
estimate its denominators we need the following property of binomial coefficients (see 
[SI lemma 4.1]): 



(22) 



if a G 



Then according to (|22|) . the numbers 

'n + a 2 - aA fn + a x - a 2 \ (a 2 + l) n+ fc 
n — / \ fc 

are integers for = 0, 1, . . . , n. This proves the lemma. 



u 2n u n 



(n + k)\ 



□ 



Remark 1. From Lemma [T] it follows that R n (a\,a 2 ,b) = —R n (a 2 ,ai,b). Note that this 
symmetry property can be obtained directly from (JIBl) . since we have Qj% 



(ai,a 2 ,b) 



X) 



Qn' a (x). The last equality follows easily from the Rodrigues formula (jHJ) and the 
fact (see [2, Ex.]) that the Rodrigues operators corresponding to the weights w± and 
w 2 in (jSJ) commute. 

Note that by the analogous arguments, the Rivoal construction ( 1T21) and therefore, 
the polynomials P n (a, x), Q n (a, x) possess a similar symmetry. Indeed, if we define the 
generalized Laguerre polynomial (see [2} §3.2]) 



X 



TV. 



{ai + l) n (a 2 + 1)„ x ( a x + n + 1, a 2 + n + 1 
^ e 2 ^ 2 



ai + 1, a 2 + 1 



— x 



10 



and corresponding to it integrals 



and 



/?„.,.(.«■)= / ^^W 'rft j = 1,2 

x — t 



S n (a 1 ,a 2 ,x) 

(23) 



then we will have A ni(11)(l2 (x) = Ai,a 2 ,aiO), ^ai.aa^) = -#n,a 2 ,ai(£)- This implies that 
the "determinantal" sequence 

-<4n,ai ,a 2 (^0 Rn,ax,a<2 (^v 
y 4n+l.ai,a2 (*^) -Rn+l.ai .02 (•^0 

= Qn(ai, a 2 , x) x a2 T(ai + 1) - -P„.(ai, a 2 , x) x ai r(a 2 + 1) 

has the following properties: 

<S n (ai, a 2 , x) = -S n (a 2 , ai, x) and P n (ai, a-2, x) = Q n (a2, a 1 ,x). 

Taking; = in 0231) we get fll2p . which yields 

Q n {a, x) = Q n (a, 0, x), P n (a, x) = Q n (0, a, x). 

It would be of interest to find a compact hypergeometric form for the polynomial 
Qn(ai, a 2 , x), if it exists (maybe as a product of hypergeometric series), for better 
understanding of hypergeometric origin of such sequences. 

The next lemma provides an alternative representation of R n (ai,a 2 ,b) as a multiple 
integral. 

Lemma 2. Let ax, a 2 , b G 1R, a 2 — ai ^ Z, a 1; a 2 > —1, 6 > 0. TTien one /ias 

1 + a 2 - ai) n (l + ai - a 2 ) n 6 2n+1 



i? n (ai, a 2 , 6) = (a 2 - a x )- 



nl 



2 



nX n+ai y n (x - l)2n+l e -te 
(X?/ + l)"+l-H»i-«a(y + l)n+l+oa-ai 35 ^ 

Proof. Substituting (jSJ) into (JIB]) we obtain 

L2n+1 /-00 i _ ai-a 2 , . 

i2 B (a 1 ,a a ,6) = ^- i =- / ^ {x n+a ^(x n+ ^(l - x) 2n+1 e~ bx ) {n) ) {n) dx. 

n! ,/ 1 — x 

Applying the following representation introduced by Rivoal (see [TUl proof of Lemma 
4]): 

1 1 - x ai " a2 f 00 du 



ai-a 2 1 - x 7 (1 + Mx) 1+a2 ~ ai (l + ii)l+ai-«a 

we get 

(ai - a 2 )6 2n+1 Z" 00 /"°° (a; n + a2 - a i(a;™ +ai (l - x) 2n + 1 e- fa )W) (n) 
R n (ax,a 2 ,b) = - / / -. — ; r— dudx. 

Now using Fubini's theorem and n-integrations by parts with respect to x we obtain 

,S {a 2 - ax) n+1 b 2n+l f°° f°° U n X n+a *- ai (x n+ai (x - l)2n+l e -fer)(n) 

Rnlax, a 2 , V) — — / / — r — — - r— dudx. 

v ' ' ; n! 2 J J (1 + wx) n+1+a2 - ai (l + u) 1+a ^ 
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Making the change of variable u — 1/ (xy) or y — 1/ (to) in the last integral we get 

Rn{ai,a 2 ,b) = / / , ■ rr— t— — 7— , rrr— — tferfj/. 



n! 2 7 7 (2/+ l) n + 1+a 2- a i(xy + l)l+ai-a 2 

Finally, integrating by parts with respect to x we obtain the desired representation and 
the lemma is proved. □ 
Our next concern will be the asymptotic behavior of the sequences q n (ai,a 2 ,b) and 
R n (a 1 ,a 2 ,b). 

Lemma 3. Let a\,a 2 ,b G R, a 2 — a% (jL Z, a±,a 2 > —1, b > 0. TTien i/ie following 
asymptotic formulas are valid: 

Rn{ax,a 2 ,b) = (2n)! ^ 1/4 _ (ai+aa)/2 (2sin(7r(a 2 - ai)) ^ c 2 + 0(n- 1/2 )), 

q n (a 1 ,a 2 ,b) = (2n)! nl/4 _ (ai+a2)/2 (c 2 + 0(n~ 1/2 )) as n 00, 

where the constant c 2 is defined in Theorem 1 . 

Proof. To obtain the asymptotic behavior of the above sequences, we proceed anal- 
ogously to the proof of Lemma 5 in [6]. Substituting (jHJ) into (TTT]) and integrating by 
parts we have 

( — '\\ n h2n+a 2 +l poo / n +a 2 -a 1 \ ( n ) 

1n(a u a 2 ,b) = 1 \ } ° ■ - / - l) 2 " +1 e^ * — - dx. 



n! r(a 2 + l)7 v ' \{x-\) n + x 

Further, repeating the arguments such as those used in [U p. 56-59] we get 

L2n+a 2 +l „2n+(ai+a 2 )/2+3/4 roc 

gn(ai,a 2 ,6)= -= / e^<»Vi(y)dl/-(l + 0(n- 1/a )), 

T(a 2 + 1) 20T J 

where 

</?!(</) = 21n(Vy + +lny - 6y, Vifo) = (Vv + l/v^)z/ (ai+a2)/2 " 1/4 - 

The last integral coincides with the corresponding integral from [6J Lemma 5] and 
likewise, by the Laplace method, we obtain 

«•<"*• 6) = filTTF^^ 6 ' 2 " /8(1 + 0( ""' /2)) 

= (2n)! n <°'+«>' 2 - 1 ' 4 e^(c 2 + Ofn-" 2 )), 

where the constant c 2 = c 2 (ai,a 2 ,6) is defined in Theorem [TJ 

Making the change of variables x = nu, y = vj 'y/rvu in the double integral of Lemma |2] 

we get 

= r(n + 1 + a 2 - aOrtn + 1 + a, - a,) ^.^^^ f» 
n! 2 r(a 2 - ai)r(l + ai - a 2 ) Jo 

where 

„. \ uv(u-l/nfe- bu . , w( ai+a2 )/ 2 (M- 1/n) 

f{u,v) = 7 — ■ — — — , </(w,u A - 



v + yjmi)(vyjmi+ 1)' + v /^) 1 + a2 - ai (w v /mZ+ l) 1+a i- a2 ' 
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Using the reflection formula 

7T 



r(a 2 - ai) • r(l + ai — 02 J 



sin(7r(a2 — ai)) 
and Stirling's formula for the gamma function 

T(n + l + a 2 -a 1 ) = n a2 ' ai n\(l + 0(n' 1/2 )) as n -4 00 

we get 

^ = sin(7r( a2 - ai )) 62w+lw3Tt+2+(ai+a2)/2 Z" 00 ^ ^ /n( ^ ^ ^ ^ + 0(n -i^ 

^ Jo 

Now applying the Laplace method for multiple integrals by the same way as in the 
proof of Lemma 5 [B] with the only difference that in our case 

n l+ a 2 

we obtain 

(r, \2n+l/4+(a 1 +a 2 )/2 

R. = 2 8 i„ W a 2 - a,)) (2n l /1+tn+a2)/2 g -*^W(l + 0( B -^) 

= (2n)\e- /}S: 'n'« + "''>' 2 - 1 "{2 S m(7,(a 2 - ai) f^ ^ c + 0(n 1 ' 2 )) 

and the lemma is proved. □ 
Now Theorem [1] follows easily from Lemmas HHSJ 

3. Proof of Theorem [2J 

It is readily seen that the construction (J7J) is a limiting case of Theorem [TJ Indeed, 
comparing Lemma [2] with Lemmas 3, 4 from [B] we get 

(24) lim - — ^h^S = Pn (a 1 ,b)-(kib-4>(a 1 + l))q n (a 1 ,b). 

02->ai 1 [CLi + 1J <2 2 — Oi 

On the other hand, by Lemma [U we have 

hm rv nT = hm [Qn{ai,a 2 ,b)— — — — - q n {a 2 , a u b) ) 

a2^ai 1 (ai + 1) a 2 — ai a2->ai a 2 — ai \ 1 (ai + l)/b ai J 

q n (a 1 ,a 2 ,b) -T(a 1 + I) q n (ai, a 2 , b) - q n (a 2 , ai, b) 
— hm — — — ■ — h hm 

a2->ai r(ai + 1) Cl 2 — 0-1 (J2^ai a 2 — <2i 

_ g w (ai,ai,&) d ^ r(a 2 + 1) ^ ^ g ra (ai,a 2 ,&) - q n (a 2 ,ai,b) 



a 2 - J >a 1 a 2 - Oi 

a2=ai * x 



r(ai + 1) da 2 V & a2_ai 
Taking into account that 

j m» + mr(., + i) 

da 2 V & a2_ai / & a2 ~ ai 
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we obtain 

(25) 

lim R ^ a ^ = gri (a 1 ,6)(^(a 1 +l)-ln6)+ lim ^{a^b) - qn {a 2 ^b) 

a2->ai T(ai + 1) a 2 — a-i a 2 ^a ± a 2 — ai 

Comparing the right-hand sides of ( 124"1) and ( 1251) we conclude that 

, ,x ,. q n (ai,a2,b) - q n (a 2 ,a 1 ,b) 
p n {a 1 ,b)= hm . 

i2->ai <2 2 — d\ 

To evaluate the above limit, we need the following expansions: 
'n + a\ — a 2 \ ( n 



k / \k 



n + a 2 — ai 
n — k 



- (H n - H n „ k )(a 2 - a x ) + 0((a 2 - aO 2 ), 
) = (k) + (^(^-^K-aO + O^-aO 2 ), 



(a 2 + l) n +k = («i + l)n+fc + (ai + l)„+fc-ffn+fe(ai)(a 2 - ai) + 0((a 2 - ai) 2 ). 
Then a straightforward verification shows that 



n / \ 2 

Pn(«i, 6) = 6 "~* \ k ) ^ + l)n+fc(^ n+*(ai) + 2#„- fc - 2H k 



and this is precisely the assertion of Theorem [2j □ 
Setting ai = 0, b = 1 yields Corollary HI 

Remark 2. Note that an analogous conclusion can be drawn for the polynomials 
P n (x),Q n (x) from (TTTj) . Comparing [9j Lemma 1] with [TOl Lemma 4] gives 



S n (x) = lim — \ ■ S n (a,x). 



Then analysis similar to that in the proof of Theorem [2] shows that Q n (x) = Q n (0, x) 
Q n (0, 0, x) and 

Q n (a,x) - P n (a,x) Q n (a,0,x) - Q n (0,a,x) 



Pr, (x) = lim = lim 



n\^J 11111 ) 

where the polynomial Q n (ai,a 2 ,x) is defined in Remark 1. 
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